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PROTON  BEAM  DEPOSITION,  HEATING,  AND 
RADIATION  FROM  AN  IRON  PLASMA 


I.  Introduction 

With  current  charged  particle  beam  technologies,  it  is  possible  to 
bombard  targets  with  intense  proton  beams.  The  interest  in  charged 
particle  beam-target  interactions  is  enhanced  by  the  expectation  that  the 
beam-target  energy  coupling  can  be  understood  ar.d  modelled,  to  first 
approximation,  in  terms  of  well  understood  classical,  linear  processes,  as 
opposed  to  the  non-linear  effects  involved  in  laser-target  coupling.1 
Thus,  much  effort  is  directed  toward  investigating  these  phenomena. 
Significant  applications  of  this  research  include  the  use  of  ion  beams  for 
x-ray  generation,  inertial  confinement  fusion,  and  weapons  lethality  and 
vulnerability  studies. 

II.  Theoretical  Model 

The  interaction  of  a  monoenergetic  one  MeV  proton  beam  with  a  planar 
aluminum  slab  of  15  um  thickness  has  been  previously  investigated.  This 
beam  was  assumed  to  consist  of  a  square-shaped  pulse  of  10  nanosecond 
duration  with  a  flux  of  102^  protons  cm-2  sec-1.  The  beam  intensity  on 
target  was  1.6x10^  W/cm2,  which  is  comparable  to  Intensities  available 
with  current  pulsed-power  generators. 

In  the  present  study,  the  interaction  of  this  same  proton  beam  with  a 
planar  8.0-um-thick  iron  slab  is  treated.  This  distance  approximately 
corresponds  zo  the  range  of  a  one  MeV  proton  in  cold  solid  density  Fe^. 

The  scheme  for  treating  the  proton  beam-iron  slab  interaction  is 
fundamentally  the  same  as  that  described  in  Ref.  2  for  treating  the  proton 
beam-Al  slab  interaction.  The  hydrodynamic  response  is  treated  via  a  one¬ 
dimensional  model  with  a  sliding-zone  version  of  flux  corrected 
transport.  A  special  gridding  algorithm  moves  zones  in  a  Lagrangian 
fashion  and  adjusts  the  mesh  in  order  to  resolve  steep  gradients  in  the 
flow.  A  single  temperature  model  is  employed,  since  equilibration  times 
are  generally  short  compared  with  hydrodynamic  timescales  in  this  problem. 
A  collisional  radiative  equilibrium  treatment  is  used  to  obtain  the  atomic 
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level  populations.  Radiation  transport  is  treated  with  a  hybrid  scheme 
which  uses  multifrequency  ray  tracing  for  the  contir.ua,  and  a  probability- 
of-escape  method  for  lines.  The  radiation  transport  is  coupled  self- 
consister.tly  into  the  ionization  dynamics.  The  proton  stopping  power 
includes  contributions  from  bound  atomic  electrons,  free  plasma  electrons, 
and  plasma  ions.  These  features  are  discussed  in  detail  in  Ref.  2. 

The  Fe  atomic  model  used  in  this  calculation  did  not  have  excited 
level  structure  in  ionization  stages  below  Fe  XVI.  Only  ground  states  were 
carried  for  stages  I  through  XV.  While  this  lack  of  level  structure  and 
the  resulting  omission  of  lines  from  the  radiation  transport  scheme  can 
modify  the  energy  balance  in  regimes  where  these  lower  stages  are  highly 
populated,  it  is  estimated  that  the  effect  is  less  than  a  factor  of  two  on 
the  energy  lost  due  to  radiation. 

A  major  difference  between  this  work  with  Fe  as  the  target  material 
and  the  previous  work  with  A1  as  the  target  material  is  in  the  treatment  of 
the  bound  electron  stopping  power.  For  the  proton-Al  +  CJ  interaction,  where 
q  is  the  ionic  charge,  detailed  stopping  power  cross  sections  are 
available;4  for  the  proton-Fe+<1  interaction,  no  such  data  are  currently 
available.  In  the  present  treatment,  stopping  power  due  to  bound  electrons 
is  calculated  from  a  method  proposed  by  Mehlho'rn.5  In  this  approach, 
stopping  power  is  calculated  from  the  Bethe  equation  (which  is  valid  at 
high  projectile  energies)  and  from  the  LSS  model  (which  is  valid  in  the  low 
projectile  energy  regime);  then  the  stopping  cross  section  is  taken  to  be 
the  lesser  of  these  two  results;  i.e., 

VE)  ‘  min  CSBethe(E)’  WE>1  ' 

Mehlhorn^  has  found  that,  by  including  shell  corrections  in  the  Bethe 
equation,  he  can  adequately  represent  cold  target  stopping  power  results. 
By  applying  this  model  to  cold  solid  density  Fe,  it  has  been  verified  that 
the  Andersen-Ziegler^  results  for  proton  stopping  power  can  be  reproduced 
very  well. 

The  Bethe  model  for  the  stopping  cross  section  is  given  by-* 
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where  A 


Here  Zg  is  the  effective  charge  of  the  projectile  ion,  Z  is  the  atomic 
number  of  the  target,  Z  is  the  average  charge  of  the  target  ior.s,  e  and  m 
are  the  electron  charge  and  rest  mass,  E  and  V  are  the  projectile  energy 
and  velocity,  3  -  V/c,  where  c  is  the  vacuum  velocity  of  light,  <I>  is 
the  average  ionization  potential  of  the  target  ions,  and  the  bi  are  least- 
squares  fit  parameters  derived  by  Andersen  and  Ziegler^  to  calculate  shell 
corrections  to  the  Bethe  stopping  power. 

Y  -  (1-S2)'1/2  . 

The  LSS  model  is  described  by  the  following  equations:^ 
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A  »  A2/A1  ,  where  A2  is  the  atomic  weight  of  the  target  ion,  and  A1  is  the 
atomic  weight  of  the  projectile  ion. 

The  effective  charge  of  the  projectile  ion  is  assumed  to  be  given  by 
the  Brown  and  Moak^  relation 

ZQ/Z1  -  1  -  1.034  exp  [-V/(zJ’69  X  Vb)] 

where  Z-,  is  the  atomic  number  of  the  projectile  ion  and  V&  is  the  Bohr 
velocity  2.188  X  10®  cm/sec. 

To  calculate  an  average  ionization  potential  for  a  target  ion  of 
charge  q,  a  scheme  proposed  by  Mehlhorn^  was  adopted; 


<1  >  -  (Z/Z)11  <1  >  , 
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where  Zp  is  the  atomic  number  and  <Ip>  is  the  average  ionization  potential 
of  a  neutral  atom  with  the  same  electronic  configuration  as  the  target  ion 
of  charge  q. 

Once  the  stopping  cross  section  is  obtained,  the  stopping  power  due  to 
bound  atomic  electrons  is  given  by 


p(x) 


Sb(E) 


where  p(x)  is  the  mass  density  of  target  ions  in  the  plasma  at  position  x, 
and  m-r  is  the  mass  of  the  target  ion. 

Stopping  power  due  to  free  plasma  electrons  and  plasma  ions  is 
calculated  from  the  same  models  as  used  in  the  previous  treatment  of  an  A1 
target.  Enhanced  stopping  power  due  to  collective  beam  effects  is 
included,  although  in  the  case  treated  here,  this  effect  is  small. 


III.  Results 


Some  results  of  the  proton  beam-Fe  slab  target  interaction  simulation 
will  now  be  given.  Figure  1  shows  the  density,  temperature,  and  average 
charge  variations  at  the  front  and  rear  surfaces  of  the  plasma  as  the 
interaction  evolves.  At  the  front  surface,  temperature  and  Z  rise  very 
quickly  and  attain  approximate  steady  state  values  within  the  first  few 
nanoseconds.  The  density  drops  rapidly  at  first,  then  drops  more  slowly  at 
later  times,  but  it  is  always  decreasing  as  the  plasma  expands.  There  is 
some  slight  cooling  and  a  noticeable  drop  in  Z  after  the  beam  is  shut  off 
at  10  nsec.  The  temperature  and  Z  rise  at  the  rear  surface  is  much  more 
gradual  and  does  not  begin  to  levei  off  until  near  the  end  of  the  beam 
pulse.  The  density  decreases  here  also  as  the  plasma  expands.  The  peaks 
in  the  rear  surface  plots  at  early  times  are  due  to  the  arrival  of  a  shock 
wave. 

Density  and  temperature  profiles  at  different  times  in  the  plasma 
evolution  are  shown  in  Fig.  2.  In  these  plots,  the  beam  is  incident  from 
the  right.  The  expansion  of  the  plasma  with  time  is  easily  discerned.  At 
early  times,  there  are  significant  gradients  near  the  rear  surface.  As 
time  goes  on,  these  variations  disappear  and  the  plasma  becomes  more 
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uniform,  especially  at  about  13  nsec,  three  nsec  after  the  beam  has  been 
turned  off.  After  three  nsec,  the  temperature  throughout  most  of  the 
plasma  remains  fairly  steady.  At  1.7  nsec,  the  beam  is  stopped  in  the 
plasma  interior  and  does  not  reach  the  backside  to  heat  it  up.  This  is 
also  true  at  later  times;  the  temperature  increase  at  the  rear  surface 
arises  from  heating  due  to  radiation  transport  from  the  interior  of  the 
plasma. 

Calculated  front  and  rear  spectra  at  these  times  are  displayed  in  Fig. 
3.  The  lack  of  level  structure  in  the  lower  ionization  stages  is 
demonstrated  by  the  paucity  of  lines  in  Fig.  3a  and  in  the  rear  spectra  in 
Figs.  3b  and  3c.  Lines  appear  in  the  spectra  when  the  plasma  gets  hot 
enough  so  that  ions  with  level  structure  have  a  significant  abundance.  The 
peaks  in  the  rear  spectra,  except  at  the  lowest  photon  energies,  arise  from 
radiation  from  the  plasma  interior  "shining  through"  opacity  "windows". 
The  dominant  contribution  to  the  front  surface  spectra  generally  comes  from 
the  outermost  zone,  although  "shine-through"  contributions  from  the 
interior  can  produce  some  peaks  in  the  curve.  As  can  be  noted  from  Figs.  1 
and  2,  the  plasma  is  nearly  uniform  at  13.15  nsec;  hence  the  front  and  rear 
spectra  are  very  similar.  The  front  spectra  becomes  more  intense  at  the 
higher  frequencies  up  to  10  nsec,  but,  as  might  be  expected,  the  13.15  nsec 
front  spectra  is  less  intense  than  at  10  nsec  due  to  plasma  cooling  after 
beam  shutoff. 

Figure  4  gives  the  energy  history  and  partitioning  of  the  beam-target 
interaction;  "dep"  is  the  energy  deposited  by  the  beam,  and  "sum"  is  the 
total  energy  remaining  in  the  plasma  in  various  forms.  After  the  beam 
shutoff  at  10  nsec,  the  total,  internal,  and  potential  energies  decrease  as 
the  plasma  cools;  only  the  kinetic  energy  continues  to  rise  as  the  plasma 
expands.  At  early  times,  the  energy  lost  by  radiation  is  insignificant, 
but  by  13  nsec,  this  represents  about  a  third  of  the  deposited  energy. 
Thus,  radiation  represents  a  significant  energy  less.  If  the  Fe  atomic 
model  had  included  level  structure  and  line  radiation  in  the  lower 
ionization  stages,  the  radiation  loss  would  be  greater,  particularly  at 
early  times.  However,  as  stated  earlier,  this  lack  of  line  structure  is 
presumed  to  have  an  effect  of  less  than  a  factor  of  two  on  the  radiation 
loss . 
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IV.  Summar; 


The  interaction  of  a  one-MeV  proton  beam  with  an  intensity  of  16 
terawatts/cm^  with  a  planar  8.0-ym-thick  iron  slab  has  been  simulated  using 
a  fully  self-consistent,  one-dimensional  hydrodynamic-ionization-radiation 
transport  code.  An  extensive  atomic-level-structure  model  of  Fe  was 
employed,  and  atomic  populations  a3  functions  of  temperature  and  density 
were  calculated  from  a  collisional-radiative-equilibrium  assumption. 
Radiation  transport  was  treated  with  a  hybrid  scheme  that  combines  a  multi¬ 
frequency  formalism  for  the  continuum  with  a  probability-of-escape  method 
for  the  lines. 

Profiles  of  the  plasma  density  and  temperature  at  various  times  in  the 
plasma  evolution,  and  corresponding  front  and  rear  emission  spectra  at 
these  times,  have  been  presented.  Front  and  rear  surface  density, 
temperature,  and  average  ionic  charge  has  been  discussed,  and  the  energy 
partitioning  among  the  various  degrees  of  freedom  has  been  shown. 

Due  to  limitations  of  computer  storage  and  other  factors,  detailed 
atomic  structure  was  only  included  for  Fe  XVI  and  higher  stages.  Since  the 
maximum  average  charge  obtained  was  about  18,  radiation  losses  are 
underestimated  throughout  the  simulation;  however,  it  is  estimated  that  the 
error  is  less  than  a  factor  of  two. 

At  these  beam  intensity  levels,  no  K-shell  emission  was  predicted,  in 
contrast  to  the  earlier  study  of  the  interaction  of  the  same  proton  beam 

p 

with  an  aluminum  target.  Many  L-shell  lines  are  present  during  the  latter 
part  of  the  simulation. 

There  has  been  much  interest  recently  in  lasing  schemes  involving 
2p^3p~2p^33  transitions  in  neonlike  systems.  One  such  transition  is  the 
2p^3p(^2.^)0''2p^33(X2,^)  1  transition  at  25^. 9A  in  Fe  XVII.  One  of  the 
reasons  for  undertaking  this  study  is  to  see  if  this  transition  can  be 
excited  by  proton  oeams.  The  atomic  model  did  not  have  detailed  jj  levels 
for  Fe  XVII;  the  3p  levels  were  combined  into  twb  groups  according  to  the 
spin  of  the  2p^(^P)  inner  shell.  Nevertheless,  it  should  be  possible  to 
make  some  estimates  of  the  feasibility  of  inducing  this  transition.  After 
about  5.5  nsec,  these  3p  levels  had  fractional  population  densities  of  “5~ 
6x1 0-^  to  ~10-^.  During  this  time,  the  ion  density  was  of  the  order  1020, 
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and  the  electron  density  of  the  order  10  .  The  temperatures  in  the  plasma 

during  these  times  averaged  about  180  eV,  with  some  hotter  zones  rising 
above  200  eV.  Under  these  conditions,  it  should  be  possible  to  see  seme 
lasing  in  this  line. 

This  treatment  is  one-dimensional;  hence,  field  effects  and  lateral 
energy  flow  have  been  omitted.  These  factors  could  alter  profile  shapes 
and  magnitudes,  and  stopping  power  in  the  plasma,  and  thus  affect  the 
results  presented  here. 
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Fig.  3  Front  and  rear  emission  spectra  at  various  times 
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'  Omaha,  Nebraska  68113 

>  ATTN:  JLKS 

!  Undersecretary  of  Defense  1  Copy 

for  RSCH  and  ENGRG 
£  Department  of  Defense 

|  Washington,  D.C.  20301 

*  ATTN:  Strategic  and  Space  Systems  (OS) 

V 

)  Deputy  Chief  of  Staff  for  RSCH  DEV  and  ACQ  1  Copy 


I  Department  of  the  Army 
|  Washington,  D.C.  20301 
.  ATTN:  DAMA-CSS-N 


Commander  i  copy  each 

Harry  Diamond  Laboratories 
Department  of  the  Army 
2800  Powder  Hill  Road 
Adelphi ,  MD  20783 
ATTN:  DELHD-N-NP 
ATTN:  DELHD-R  J.  Rosado 
ATTN:  DELHD-TA-L  (Tech.  Lib.) 

U.S.  Army  Missile  Command 

Redstone  Scientific  Information  Center  3  Copies 

Attn:  DRSMI-RPRD  (Documents) 

Redstone  Arsenal,  Alabama  35809 

Commander  1  copy 

U.S.  Army  Missile  Command 
Redstone  Arsenal,  Alabama  35898 
ATTN:  DRCPM-PE-EA 

Commander 

U.S.  Army  Nuclear  and  Chemical  Agency 

7500  Backlick  Road  1  copy 

Building  2073 
Springfield,  VA  22150 
ATTN:  Library 

Commander  1  Copy 

Naval  Intelligence  Support  Center 
U301  Suitland  Road,  Bldg.  5 
Washington,  D.C.  20390 
ATTN:  NISC-U5 

Commander  1  Copy 

Naval  Weapons  Center 
China  Lake,  California  93555 
ATTN:  Code  233  (Tech.  Lib.) 


Officer  in  Charge  1  Copy  each 

White  Oak  Laboratory 
Naval  Surface  Weapons  Center 
Silver  Spring,  Md.  20910 
ATTN:  Code  R40 
ATTN:  Code  F31 

Air  Force  Weapons  Laboratory  1  Copy  each 

Kirtland  AFB,  New  Mexico  87117 


ATTN: 

SUL 

ATTN: 

CA 

ATTN: 

APL 

ATTN: 

Lt.  Col  Generosa 

1  Copy 


Deputy  Chief  of  Staff 
Research,  Development  and  Accounting 
Department  of  the  Air  Force 
Washington,  D.  C.  20330 
ATTN:  AFROQSM 

Commander  1  Copy 

U.S.  Army  Test  and  Evaluation  Command 
Aberdeen  Proving  Ground,  MD  21005 
ATTN:  DRSTE-EL 

AVCO  Research  and  Systems  Group  1  Copy 

201  Lowell  Street 
Wiimi.nton,  MA  0’ 387 

ATTN:  Library  A330 

3 DM  Corporation  1  Copy 

7915  Jones  Branch  Drive 
McLean,  Virginia  22101 

ATTN:  Corporate  Library 

Berkeley  Research  Associates  1  Copy 

P.0.  Box  983 
Berkeley,  CA  9^701 

ATTN:  Dr.  Joseph  Workman 

Berkeley  Research  Associates  1  Copy  each 

P.0.  3ox  852 
5532  Hempstead  Way 
Springfield,  VA  22151 

ATTN:  Dr.  Joseph  Orens 

3oeing  Company  1  Copy 

P.  0.  Box  3707 
Seattle,  WA  98 1 3 u 

ATTN:  Aerospace  Library 

The  Dikewood  Corporation  1  Copy 

1613  University  Bldv.,  N.E. 

Albuquerque,  New  Mexico  87110 
ATTN:  L.  Wayne  Davis 

General  Electric  Company  1  Copy 

Space  Division 

J  Valley  Forge  Space  Center 

P.  0.  Box  8555 
Philadelphia,  PA  19101 
ATTN:  J.  Peden 
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General  Electric  Company  -  Tempo 
Center  for  Advanced  Studies 
816  State  Street 
P.0.  Drawer  QQ 
Santa  Barbara,  CA  93102 
ATTN:  DASIAC 

Institute  for  Defense  Analyses 
1801  N.  Beauregard  St. 
Alexandria,  VA  22311 

ATTN:  Classified  Library 

IRT  Corporation 
P.3.  Box  81087 
San  Diego,  CA  92138 
ATTN:  R.  Mertz 


1  Copy 


1  Copy 


1  Copy 


JAYCOR 

11011  Forreyane  Rd. 
P.0.  Box  851 5^ 

San  Diego,  CA  92138 
ATTN :  E .  Wenaas 
F.  Felbar 


1  Copy 


JAYCOR 

205  S.  Whiting  Street,  Suite  500 
Alexandria,  VA  22304 
ATTN:  R.  Sullivan 


1  Copy 


KAMAN  Sciences  Corp. 
P.  0.  Box  7463 
Colorado  Springs,  CO 
ATTN:  Library 


1  copy  each 


80933 


Lawrence  Livermore  National  Laboratory 
University  of  California 
P.0.  Box  808 

Livermore,  California  94550 
Attn:  DOC  CDN  for  L-153 
Attn:  DOC  CDN  for  L-47  L.  Wouters 
Attn:  DOC  CDN  for  Tech.  Infor.  Dept.  Lib. 

Lockheed  Missiles  and  Space  Co.,  Inc. 

P.  0.  Box  504 
Sunnyvale,  CA  94086 
Attn:  S.  Taimlty 
Attn:  J.D.  Weisner 

Lockheed  Missiles  and  Space  Co.,  Inc. 

3251  Hanover  Street 
Palo  Alto,  CA  94304 
Attn:  J.  Perez 


1  copy  each 


1  copy  each 


1  Copy 
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1  Copy  each 


Maxwell  Laboratory,  Inc. 

9244  3alboa  Avenue 
San  Diego,  CA  92123 
ATTN:  A.  Kolb 
ATTN:  M.  Montgomery 
ATTN :  J .  Shannon 

McDonnell  Douglas  Corp. 

5301  Solsa  Avenue 
Huntington  3each,  CA  9264? 

ATTN:  S.  Schneider 

Mission  Research  Corp. 

P.  0.  Drawer  719 
Santa  3aroara,  CA  93102 
ATTN:  C.  Longmire 
ATTN:  W.  Hart 

Mission  Research  Corp. -San  Diego 
5434  Ruffin  Rd. 

San  Diego,  California  92123 
ATTN:  Victor  J.  Van  Lint 

Northrop  Corporation 

Northrop  Research  and  Technology  Center 
1  Research  Park 

Palos  Verdes  Peninsula,  CA  90274 
ATTN:  Library 

Physics  International  Company 
2700  Merced  Street 
San  Leandro,  CA  94577 
ATTN:  M.  Krishnan 
ATTN:  C.  Gilman 
ATTN:  S.  Wong 

R  and  D  Associates 
P.0.  Box  9695 
Marina  Del  Rey,  CA  90291 
ATTN:  W.  Graham,  Jr. 

ATTN:  P.  Haas 

Sandia  National  Laboratories 
P.0.  Box  5800 

Albuquerque,  New  Mexico  87115 


ATTN: 

Doc  Con  For  3141 

ATTN: 

D.  McDaniel 

ATTN: 

P.  VanDevender 

ATTN: 

K.  Matzen,  Code  4247 

Science  Applications,  Inc. 

P. 

0.  Box 

2351 

La 

Jolla, 

CA  92038 

ATTN: 

R.  Beyster 

1  Copy 


1  Copy  eac.n 


1  Copy 


1  Copy 


1  Copy  each 


1  Copy  each 


1  copy  each 


1  copy 
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1  copy 


Spectra  Technol,  Inc., 

Attn:  Alan  Hoffman 
2755  Northup  Way 
Bellevue,  Washington  98004 

Spire  Corporation  1  copy 

P.  0.  Box  D 
Bedford,  MA  01730 
ATTN:  R.  Little 

SRI  International  1  copy 

333  Ravenswood  Avenue 
Menlo  ParK,  CA  94025 
ATTN:  S.  Dairiki 

S-CUBED  1  copy 

P.  0.  Box  1620 
La  Jolla,  CA  92038 
ATTN:  A.  Wilson 


Director 

Strategic  Defense  Initiative  Organization 
1717  H  Street 
Pentagon  20301-7100 

ATTN:  DE  Lt .  Col  Richard  Guliickson/DEO 
1ST  Dr.  Dwight  Dus ton 
ATTN:  1ST  Dr.  J.  Ionson 

Texas  Tech  University 
P.0.  Box  5404 
North  College  Station 
Lubbock,  TX  79417 
ATTN:  T.  Simpson 

TRW  Defense  and  Space  Systems  Group 

One  Space  Park 

Redondo  Beach,  CA  90278 

ATTN:  Technical  Information  Center 


University  of  Buffalo 

Attn:  Professor  W.  J.  Sarjeant 

Dept,  of  Electrical  Engineering 

High  Voltage  and  Power 

312  Bonner  Hall 

Buffalo,  New  York  14260 


Naval  Research  Laboratory 
Plasma  Radiation  Branch 


Washington,  D.C 
Code  4720 
Code  4700 
Code  2628 


20375 

50 

Copies 

26 

Copies 

22 

Copies 

1  copy  each 


1  copy 


1  Copy 


1  Copy 
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